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Abstract 

Muon-induced neutrons constitute a prominent background component in a number of low count rate 
experiments, namely direct searches for Dark Matter. In this work we describe a neutron detector to measure 
this background in an underground laboratory, the Laboratoire Souterrain de Modane. The system is based 
on 1 m'^ of Gd-loaded scintillator and it is linked with the muon veto of the EDELWEISS-II experiment for 
coincident muon detection. The system was installed in autumn 2008 and passed since then a number of 
commissioning tests proving its full functionality. The data-taking is continuously ongoing and a count rate 
of the order of 1 muon-induced neutron per day has been achieved. 
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1 Introduction 

A common challenge in rare-event searches is an ac- 
curate and precise understanding of all background 
components. A typical approach is based on a com- 
bination of Monte-Carlo simulations and dedicated 
measurements. In this work we present a detection 
system to study the production of muon-induced neu- 
trons in an underground laboratory. The main goal is 
to perform a quantitative analysis of this background 
component in the specific environment of Laboratoire 
Souterrain de Modane (LSM) within the context of 
EDELWEISS-II m [2] and EURECA [3] direct dark 
matter (DM) search experiments. Such experiments 
probe nature for weakly interacting massive particles 
(WIMP), as the most promising DM candidate, by 
means of nuclear recoils originating from a WIMP- 
nucleon elastic scattering. The expected interaction 
rate of WIMPs is extremely small, i.e. below 0.01 
event per day and per kilogram of the target mate- 
rial, thus rising up the importance of the background 
knowledge. The recent developments in discrimina- 
tion techniques of the DM detectors [3] allow a very 
powerful rejection of the electromagnetic background 
while neutrons reaching a detector with energies of 
the order of 1 MeV, in the case of EDELWEISS-II, 
can cause a signal potentially indistinguishable from 
the WIMP-nucleon scattering. The neutrons arise 
in an underground laboratory either from natural 
radioactivity having initial energies below 10 MeV, 
which can be effectively shielded, or being produced 
in a muon-induced cascade with energies up to several 
hundred MeV and able to travel far from a produc- 
tion vertex. In this paper we concentrate on studying 
the latter ones. The effect of the muon-induced neu- 
tron background can be reduced considerably by tag- 
ging the original muons with a muon veto around an 
experiment. However, the neutrons produced in the 
rock by muons passing outside of the setup and hence 
missed by the veto, can still reach the DM detectors. 
These neutrons can also lead to an in-situ produc- 
tion of secondary neutrons via spallation, especially 
in a lead (Pb) material commonly used around the 
DM detectors for an efficient gamma attenuation. In 



addition, a fraction of muons can be missed due to a 
non-ideal muon veto efficiency and result in a neutron 
production inside the shielding, again predominantly 
in the Pb material. 

In general, to measure the muon-induced neutron 
background is a challenging task due to its low in- 
tensity and its dependence on the actual experimen- 
tal set-up. For instance, one could look for coinci- 
dences between the muon veto and the DM detec- 
tors but as the current data-taking of EDELWEISS- 
II shows, only few muon-induced events are identi- 
fied in an exposure of several hundred kg. days [21 E]- 
To perform corresponding Monte-Carlo simulations 
is not a simple issue either, since many parameters 
have to be taken into account, e.g. the muon en- 
ergy spectrum and the fiux in an underground labo- 
ratory, production of neutrons by hadronic and elec- 
tromagnetic muon-induced showers, material depen- 
dency of this production, three dimensional topology 
of interactions. As reviewed in [JllEj, measurements 
of neutron production rates in liquid scintillator at 
different laboratory depths can be well described by 
most frequently used packages, GEANT4 |U] and 
FLUKA [lOj, which are in agreement with each other 
within <30% at the muon energies specific for an un- 
derground laboratory. However, the calculated neu- 
tron yield for Pb and an ordinary rock can vary up 
to a factor of 2 between the packages, being lower 
in GEANT4 [7]. Experimental neutron production 
rates in Pb measured in [HI |T2] and scaled to the 
same mean muon energy appear to be inconsistent 
within a factor of 3 among each other and may sug- 
gest that both simulation packages underestimate the 
neutron yield in heavy targets [3 [T3] . On the other 
hand, recent studies performed in [T^ with an ex- 
isting experimental setup though not fully optimized 
to the specific needs of muon-induced neutron detec- 
tion result in a measured neutron yield in Pb being 
1.8 (3.6) times smaller than predicted in GEANT4 
(FLUKA) . Overall, these measurements give a spread 
of neutron production rate in Pb of nearly one order 
of magnitude which can affect the sensitivity of DM 
experiments by the same factor. 

In order to measure the muon-induced neutron 
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background and to develop the corresponding back- 
ground model optimized for the LSM laboratory 
hosting the EDELWEISS-II experiment and the fu- 
ture one, EURECA, we have designed and installed 
a dedicated neutron counter. Having the neutron de- 
tector in-location of present and planned DM set- 
ups allows to avoid several uncertainties in the sim- 
ulations, e.g. arising from calculations of the muon 
energy spectrum and the muon flux in the labora- 
tory. The goal is to investigate the muon-induced 
neutron yield in Pb either with or without tagging 
the incident muon to an accuracy of better than 
25% with equal statistical and systematic uncer- 
tanties expected. One should also note here, that 
the EDELWEISS-II strategy on neutron background 
investigation is complemented by continuous moni- 
toring of the thermal neutron flux in the vicinity of 
the DM experiment by means of ^He detectors |15j . 
Another new experiment to study muon-induced neu- 
tron flux deep underground was proposed by other 
authors in [16]. The investigation of the neutron 
background is, in fact, not limited to the field of 
dark matter and is relevant, for instance, for experi- 
ments looking for neutrinoless double beta decay, e.g. 
NEMO-3 dZ] and SuperNemo [IH] with their location 
in LSM too. Furthermore, the results of our studies 
of muon-induced neutrons are as well applicable to 
other underground laboratories. 



The structure of the paper is as follows: in Sect.[2j 
a general overview of muon interactions and produc- 
tion of secondary neutrons is given, emphasizing the 
role for DM search. Sect. [3] describes the details of 
the neutron detector installed in LSM, including the 
data acquisition system, LED-based control of the 
liquid scintillator and a safety monitoring. Results of 
commissioning runs demonstrating the performance 
of the system are discussed in Sect.|4j Sect.[5]presents 
preliminary data from first physics runs as proof-of- 
principle. In Sect. [6j we conclude with a summary 
and an outlook of the anticipated data taking neces- 
sary to get full physics results. 



2 Muon induced neutrons 

Atmospheric muons are mainly produced in decays 
of charged pions and kaons created in an interaction 
of primary cosmic rays with the atoms of the Earth's 
atmosphere. The energy spectrum of muons extends 
above 1 TeV while the mean energy is about 4 GeV 
at the ground |19j . The overall angular distribution 
of muons at the sea level is oc {9 being the 

zenith angle) , at lower energy it becomes increasingly 
steep, while at higher energy it flattens, approach- 
ing a l/cos0 distribution for E"^ ^ 115 GeV and 
9 < 70° |19j . The total muon energy spectrum in 
the underground lab can be calculated rather well in 
the assumption of flat overburden [Sj [20] and it be- 
comes harder, e.g. {Ej;^^'^) = 260 - 290 GeV. To 
get a corresponding angular distribution for an un- 
derground site, a mountain profile has to be taken 
into account and thus it becomes laboratory specific. 
When passing through the rock, the muons lose their 
energy in four main processes: ionization (including 
the production of high-energy (5-electrons) , produc- 
tion of e~e^-pairs, bremsstrahlung and deep inelastic 
scattering (DIS). These processes naturally also lead 
to a development of electromagnetic and hadronic 
showers. Therefore, under the muon-induced back- 
ground one has to consider, in general, all particles 
created both directly by original muons and in sec- 
ondary cascades. As it was already explained in the 
introduction, neutrons constitute the most danger- 
ous background component for direct DM searches 
among other muon-induced particles. GEANT4 sim- 
ulations show (Fig. [T]) that most of the neutrons are 
created in the secondary cascades, predominantly by 
photonuclear interaction, pion and neutron inelastic 
scattering, pion absorption at rest and proton inelas- 
tic scattering [71 [20j . The contribution of direct muon 
nuclear interaction being the second most important 
one at low muon energies decreases with higher ener- 
gies because of larger electromagnetic and hadronic 
cascades. Neutrons produced in electromagnetic in- 
teractions have typically low energy and can be ef- 
fectively absorbed. In contrast, the hadronic showers 
and especially high energy muon DIS reactions lead 
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to a production of energetic neutrons which can travel 
far away from the creation point, reach the DM set- 
up, penetrate the entire shielding and cause nuclear 
recoils in the DM detectors either directly or by a 
neutron spallation in the high-Z material surround- 
ing the experimental setup. It also has to be noted 
that neutron production generally increases with the 
muon energy and the average atomic weight of the 
material, e.g. in Pb there will be ~100 times more 
neutrons created compared to CnH2n and ~10 times 
more compared to rock or concrete of the same thick- 
ness El El HO] . To model the effect of the muon- 
induced neutrons in application to DM searches, it is 
thus essential to take into account all the processes 
aforementioned, track each single particle through 
different materials down to an energy deposit in DM 
detectors on an event-by-event base. Monte-Carlo 
simulations based on the GEANT4 toolkit (version 
S.l.pOl) were developed in [20J and included all these 
aspects, profile of the mountain, full geometry of 
the EDELWEISS-II experiment and topology of an 
event. The EDELWEISS-II set-up, described else- 
where PP, hosts in a single cryostat an array of 
Ge bolometers used as DM detectors. The cryo- 
stat is surrounded by 20 cm lead shield, followed 
by 50 cm polyethylene and a modular muon veto 
O EI]. As the simulations show [20j, most of neu- 
trons leading to energy deposit in Ge bolometers are 
due to secondary production in the Pb shield near 
the cryostat. The rate of this background for the 
optimized EDELWEISS-II set-up can be as low as 
r ~ lO'Vkg/day once the anti-coincidence with the 
muon veto and with other bolometers is required but 
crucially depends on the shielding concept and the 
muon veto efficiency. As follows from these simu- 
lations and the ones described in |8j, these muon- 
induced neutrons are not a matter of concern for the 
current phase of EDELWEISS-II aiming for the sen- 
sitivity of 10~^ pb for the spin-independent WIMP- 
nucleon scattering cross-section but may become so 
for the next phase and especially for the next genera- 
tion DM experiments with the goal of lO"^*' pb, such 
as EURECA. With the neutron counter described 
here we intend to gather statistically significant data 
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Figure 1: Relative contribution of individual processes to the 
total neutron production yield as a function of the muon energy 
in CnH2n, GEANT4 simulations. Figure is adopted from [20] . 
the lines are drawn to guide the eye. 



for validation and further development of our back- 
ground model [20] . 



3 Experimental set-up 

3.1 Measurement principle 

The LSM underground lab, where the designed neu- 
tron detector is installed, is situated in the Frejus 
road tunnel in the French-Italian Alps and has an 
overburden shielding of 4850 m.w.e. This reduces 
the muon flux down to about 5 /m^/day [22], i.e. 
by more than a factor of 10^ comparing to sea level. 
Such a low rate naturally leads to a low flux of muon- 
induced neutrons and requires a relatively large size 
of the instrument. To measure the neutron yield orig- 
inated from well-identified muons interacting in Pb, 
the developed neutron counter (NC) has a layer of Pb 
and its read-out system is coupled to the muon veto 
system of the EDELWEISS-II experiment for a good 
tagging of incident muons. In order to achieve an 
efficient detection of neutrons, a Gadolinium-loaded 
liquid scintillator is used. The primary trigger in this 
case is associated with the muon itself and follow- 
ing triggers correspond to secondary neutrons pre- 
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dominantly created in Pb and captured thereafter on 
Gadolinium (Gd). A measurement of the neutron 
production rate without demanding a primary muon 
tag can only be done by requiring a multiple trigger 
occurrence in the NC within a predefined acquisition 
window {tuAQ ~ 60 /is). Together with an appropri- 
ate energy threshold this allows to effectively reject 
accidental background. For these events the primary 
trigger arises from any particle produced in muon cas- 
cades. For instance, high energy gamma's can first 
lead to Compton scattering in the NC scintillator vol- 
ume and only then produce secondary neutrons in 
the Pb layer via photonuclear reaction. High energy 
neutrons could create a prompt recoil proton at first 
and then generate secondary neutrons being charac- 
terized as a multiple neutron capture event. In this 
mode we study the overall neutron production rate 
without identifying the primary particle. The events 
with multiplicity 2, i.e. a primary trigger and one 
secondary, may in part contain a contribution of cor- 
related backgrounds, e.g. the fast neutrons originat- 
ing from natural radioactivity or /3-decays of ^^^Bi 
followed by an a-decay of ^^^Po: The flux of the am- 
bient neutrons is 10~^ neutrons/cm^/s above 1 MeV 
[2H], the energy spectrum extends to 10 MeV but is 
strongly suppressed above 4 MeV [23l [2l] • Neverthe- 
less, these neutrons can still lead to recoil protons 
above the NC threshold and then be captured in the 
scintillator. -^^^Bi comes from the Uranium-chain, 
undergoes the /3-decay with the endpoint energy of 
3.27 MeV to ^^^Po which then a-decays with a half- 
life of 164.3 fis and E„ = 6.6, 6.9, 7.68 MeV. Thus 
both /3- and a-particles can deposit energy above 
the NC threshold within the acquisition window of 
toAQ ~ 60 fis. This implies that events with one 
secondary hit need special attention in the analysis 
while the effect of both background components is 
negligible for higher multiplicities in our studies. 

3.2 General overview 

The core of the detector consists of 1 m^ of liq- 
uid scintillator (50x100x200 cm^) loaded with 0.2% 
of Gd by weight, St. Gobain Bicron BC525. The 



neutron capture on Gd has the largest cross-section 
(o-^ = 2.54 • 10^ b for ^^^Gd, 15.65% of natural abun- 
dance |25) ) compared to other stable elements and re- 
sults in several 7's with a summed energy of 8 MeV. 
Adding 0.2% of Gd in organic scintillator leads to a 
capture time, Tcapture-, of about 17 /is while a typical 
time needed for neutrons to thermalize is defined by 
the proton density in the liquid and is ~10 (see 
Sect. 4.3). Underneath the detector, a 10-cm thick 



layer of lead bricks is put for studies of the neutron 
yield. The thickness of this layer was chosen by con- 
sidering the process of secondary neutron production, 
the probability for these neutrons to reach the scintil- 
lator and the mechanical stability of the instrument. 
The geometry of the NC scintillator volume was op- 
timized for the efficiency of neutron detection with 
specially adapted Monte-Carlo simulations based on 
GEANT4 (version 8.1.p01) and taking into account 
space and security restrictions of the underground 
laboratory. In the simulations it was found that most 
of the secondary neutrons are captured after ~20 cm 
of penetration into the scintillator. The registration 
efficiency of gammas produced by Gd(n,7) saturates 
once a 50-cm thick system is realized (Fig. [2]), thus 
defining the thickness of the instrument. However, 
this efficiency is significantly less than 100% because 
of a partial escape of neutron-capture gammas as 
a consequence of a limited neutron penetration into 
the scintillator, as well as due to an applied energy 
threshold. As an example, for En = 10 MeV and 
a threshold of 3 MeV it reaches a value of ~33%. 
For higher neutron energies, the inelastic scattering 
process (n, 2n) starts to play a role, thus increasing 
the detection efficiency. Due to the partial 7 escape 
it is obvious that, instead of the 7-peak at 8 MeV 
sum energy originating from the neutron capture on 
Gd, a continuous spectrum extending to 8 MeV is ex- 
pected. The detector volume is divided in three com- 
partments: a central one containg the 1 m^ of scintil- 
lator, and two smaller ones on either side to hold the 
photomultipliers and filled with highly transparent 
mineral oil (paraffin) (Fig. [3]). There are in total 16 
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PMTs (8 on each side) of Hamamatsu R5912 type|^ 
8-inch diameter, the gain and the read-out electron- 
ics of which are optimized for the neutron capture 
events. In addition, 6 smaller PMTs (3 on each side) 
of Philips VALVO XP 2262 type|fl 2-inch diameter, 
are installed in between of these 8-inch photomult- 
pliers. Since a typical muon signal is well above the 
dynamic range of the PMT's used for the neutron 
registration, the 2-inch PMT's are operated at a gain 
adapted to the typical muon energy deposit of the 
order of 100 MeV. The entire container holding all 
liquids is made of a transparent plexiglass. A col- 
lection of scintillation light thus benefits from the 
total reflection at the outside surfaces. In order to 
avoid a pressure build-up possible due to a tempera- 
ture variation, an expansion bag filled with an argon 
gas is connected to the scintillator volume through a 
siphon system. The plexiglass chamber is wrapped 
with an aluminum foil to further improve the light 
collection. It is then placed in an aluminum vessel as 
a safety container. Finally, the system is surrounded 
by iron plates (excluding top plane) to reflect a frac- 
tion of neutrons back to the scintillator. To allow a 
coincidence measurement of the incident muon, the 
NC is coupled electronically with the existing muon 
veto of the EDELWEISS-II experiment. This veto 
system (//Veto) is constructed from 42 plastic scintil- 
lator modules (BC-412 type from St.Gobain) with a 
total of 100 m^ of surface. It is further extended by an 
additional module (5x65x315 cm^) placed on top of 
the neutron detector and the 2-inch PMTs of the NC, 
described above. The complete set-up is positioned 
right near the western wall of the EDELWEISS muon 
veto system. 

3.3 Data acquisition 



As was described in Sect. 3.1, two event signatures 
for muon-induced neutrons are considered: i) A high 
energy muon triggers either the muon system only or 
the muon veto and the neutron counter together. A 
delayed signal due to a neutron capture is then looked 
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Figure 2: Neutron detection efficiency, i.e. number of regis- 
tered neutrons vs. number of incident neutrons, as a function 
of scintillator thickness (GEANT4 simulation). Neutron is con- 
sidered to be registered when it is captured and the following 
7's deposit in the scintillator a sum energy above 3 MeV in this 
case. Neutrons start isotropically with different energies (in- 
dicated on the plot) and equally distributed over a horizontal 
bottom surface of the detector. 



for in the NC; ii) Multiple trigger occurrence hap- 
pens in the NC detector within a defined time window 
(tDAQ ~ 60 fis). The developed data acquisition sys- 
tem (DAQ) accounts for both of these possibilities 
and allows to configure two different measurement 
subsystems separately optimized for cases i) and ii). 
It is based on VME electronics, integrated into the 
/uVeto DAQ and uses specially built software in a 
Linux environment. First, in order to enhance the 
light collection two neighboring neutron PMT's (de- 
noted A and B, Figjsf^rz^/ii) and Fig.Q are combined 
into one group with individual high voltage supply. 
Thus, already summed signals are used for an event 
selection and one measurement subsystem consists of 
4 PMT groups, 2 on each side. The signals of each 
group are split and sent to a discriminator, CAEN 
V895, as well as with a delay to a flash analog-to- 
digital converter (ADC), CAEN VX1720, 250 MS/s. 
Once a threshold is exceeded, logic ECL signals gen- 
erated by the discriminator are delivered to a in- 
home developed Logic Unit (LU), based on field- 
programmable gate array (FPGA). This LU takes 
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Figure 3; Left: general layout of the neutron counter detector, positions of LED's are indicated by stars. Right: a view from 
inside of the installed but not yet filled scintillator chamber, 8 of the 8-inch and 3 of the 2-inch PMTs are seen. Numbers 
indicate different PMT groups consisting of two 8-inch PMTs (A,B) each. 



care of the final acceptance of the event, where the 
event is considered to happen when at least two op- 
posite PMT groups related to the same measurement 
subsystem get a light signal above threshold. In the 
LU module this is realized as follows. Signals from 
two neighboring PMT groups primarily pass through 
an OR-logic. Two Oi?-units having as an input the 
signals of opposite PMT groups are then wired to an 
^iVZ?-logic, which are thus 2 in total and define our 
two measurement subsystems. A coincidence on any 
of these two AND-units gives a trigger to the flash 
ADC to convert an analog signal. After a defined 
delay of Idaq ~ 60 /xs, the LU generates a common 
stop signal, i.e. ends the acquisition cycle, actuates 
the DAQ software and withdraws a time stamp of 
the triggering event from an external clock (10 /is 
precision). The latter allows to get a direct corre- 
spondence between NC and fiVeto events. Each of 
the described above OR- and ^iVZ)-outputs is con- 
nected to a time-to-digital converter (TDC), LeCroy 
1176, which allows to have a time information of in- 
dividual hits within the event with the resolution of 
1 ns. The flash ADC allows to reconstruct an energy 
signal of every individual hit within t_DylQ without 
creating a dead time. In addition, the LU card is 
triggered also in case of a muon-like event registered 



by the /uVeto. While the DAQ software records an 
acquired event, the VME bus is locked to avoid dis- 
tortion of the data. Once the event is processed, the 
VME bus is reset and the system is ready for a next 
event. The dead time to treat an event was mea- 
sured with a pulse generator and was found to be 
6.76 lb 0.08 ms, while the typical raw count rate of 
the NC is 10 Hz. 

3.4 Monitoring systems 

When dealing with liquid scintillators, especially Gd- 
loaded ones, it is important to control their aging as 
it can influence the registration efficiency with time. 
For this purpose the NC is equipped with a dedicated 
system of light-emission-diodes (LED) of 425 nm 
wavelength, corresponding to the maximum emission 
of the scintillator. There are 8 LED's in total: seven 
located on the same flank of the NC and one, marked 
"0" , in the center of the top plain (Fig.[3](Ze/i)). Their 
control is performed from the DAQ computer via a 
specifically designed VME-based module. It is possi- 
ble to adjust the voltage of each LED and an individ- 
ual time of firing, while the duration of a light pulse 
is fixed by electronic components resulting in a PMT 
signal of about 10 ns (FWHM). The corresponding 
electronics is based on the LED calibration system 
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Figure 4: General layout of the VME-based data acquisition, see Sect. 13.3] for details 



designed for the Pierre Auger Observatory |26j. The 
control module also generates at the time of LED ig- 
nition a logic pulse, or LED flag, which is then given 
as an input to the TDC of the DAQ (Fig. This 
flag allows to clean out the data from the LED en- 
tries while searching for the muon-induced events in 
an offline analysis. The LED's are typically ignited 
three times per day with a frequency of 1 Hz in a 
predefined cycle. 

Another issue of the liquid scintillator used is its 
pseudocumene base, which requires a careful han- 
dling. A safety system was thus designed and in- 
stalled. It contains several vapour detectors, two liq- 
uid leak sensors, temperature measurement outside 
and inside of the system as well as monitoring of 
the siphon level. All parameters are accessible on- 
line through a web-based interface. More details on 
this slow control monitoring can be found in |27j . 



4 Commissioning 

The NC detector was installed in LSM in autumn 
2008. After a commissioning period, data-taking 
started in early 2009. 



4.1 Control of the scintillator trans- 
parency 

Data acquired each day with the LED's are ana- 
lyzed to monitor the light properties of the scintil- 
lator in the following way: one of the central LED's, 
e.g. LEDq, is used first to normalize the gain of the 
PMT's and to account for any temporal fluctuations 
of PMT's and electronics. Then for the measurement 
performed with a non-central LEDj, a ratio of ADC 
signals of "near" PMT's, i.e. those located closer to 

'far" PMT's is calculated: 



the chosen LEDj , and 



ADC^ 



near PMT' 



's I ADC' 



(0) 



near PMT's 



^^^far PMT's I ^-^^far PMT'. 



1..7 



Any changes of r]i over time can thus be attributed 
only to changes in light transport, predominately in 
the transparency of the scintillator. The parame- 
ters rji are monitored on an everyday basis and 4- 
months trends measured with two non-central but 
symmetrically-positioned LED's, LED2 and LED4 
(Fig. [sj^/e/i)), show consistent relative changes of the 
?7-ratio (Fig. [5]). Analysis of the combined data yields 
a degradation of 15.5ib3.1 % per year as the weighted 
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Figure 5: Long term behaviour of the scintillator transparency: 
4-months trends of the rj ratio for two non-central LED's; 



LED2(top) and LED 4,{bottom) (see Sect. 4.1 for details). Both 



data sets are fit independently with a linear function. 

average. A degrading light signal can be compen- 
sated by increasing the PMT sensitivity accordingly, 
e.g. by raising the high voltage of the PMTs. In such 
case this degradation is acceptable for an expected 
run period of the NC of 2-3 years. The fluctuations of 
the 77-ratio seen in Fig. [5] are typically correlated with 
variations of the scintillator temperature. The same 
?7-ratio is also reconstructed in GEANT4 simulations 
of the NC with light propagation included. The alu- 
minum foil is described with the groundfrontpainted 
model for a diffuse reflection. A value of 0.9 for Al 
reflectivity is taken, optical properties of the scintil- 
lator and the plexiglass are considered. A difference 
between the absolute T/-values of two measurements 
as seen in Fig. [5] can be reproduced in the simula- 
tions and can be explained by a misalignement of 
the LED's and the directional property of the LED 
light. At present, this simulation agrees with the 
measurement within 20% and further improvement 
of the light tracking model is ongoing. 

4.2 Identification of muon events 

Since the described detector is accompanied by two 
complementary systems to identify in particular a 



muon crossing the NC (/iVeto plastic scintillator 
module on top and the muon PMT's inside of the 
NC), their performance has to be checked as well. 
The average energy of muons passing the mountain 
overburden is {E^^^') = 260 - 290 GeV, thus one 
expects energy losses in a 50 cm thick detector fol- 
lowing a Landau distribution. By selecting events 
registered by both muon PMT's and the top /uVeto 
module and reconstructing their energy deposit in 
the volume of the neutron detector, one finds a good 
agreement with a Landau expectation (Fig. [6]). This 
spectrum is rather broad due to a superposition of 
various track lengths ranging from more than 50 cm 
to very short ones for grazing muons in the NC. 
The count rate of events that trigger both systems 
is = 5.1 lb 0.2 /day, from which one can deduce 
the muon flux in LSM using the geometrical area 
and the solid angle of the system. The geometri- 
cal overlap is S=1.3 m^, the solid angle is less than 
2-K due to a 15 cm gap between the /iVeto module 
and the liquid scintillator level in the NC due to the 
mechanical support. In addition, a fraction of muons 
is missed due to short track lengths of grazing muons 
leading to energy deposits below the threshold (see 
the low energy cutoff in Fig. [6]). All together, this re- 
sults in an effective reduction of the measured muon 
flux estimated to be as large as 30%. Correcting for 
this, we get an estimated muon flux of 5 /m^/day 
which is in agreement with an earlier measurement 
of ^>/, = 4.98 ± 0.09 /mVday [22]. 

4.3 Verification of neutron detection 

An energy threshold of the NC has to be set such 
that most of an intrinsic background is excluded but 
the registration efficiency for 7's after the neutron 
capture is still high enough. Hence, this value was 
chosen to be around 2.5 MeV in the center of the de- 
tector. For the measurement of the neutron yield due 
to untagged muons the energy threshold is set to an 
even higher value of 3.5 MeV to further suppress ac- 
cidental and correlated background. In order to put 
a proper value of the electronic threshold and to tune 
the gain of neutron PMT's, a ^'^Co source was used 
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Figure 6: Energy deposition measured with the muon PMT's 
of NC for the events in coincidence with the /xVeto module 
on top of the NC, live time is 177.9 days. Last bin contains 
events with energy deposition of E>440 MeV (overflow bin). 
Landau distribution is fitted to the data and calibrated with 
MC simulations of muons passing through the NC. 
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Figure 7: Energy spectra acquired with AmBe source (spectra 
of the first and second ADC hits are shown) [solid lines). Both 
ADC hits of AmBe measurement are compared with corre- 
sponding GEANT4 simulations {crosses, dashed line){see text 
for details). 



first. Two gammas of 1173.2 keV and 1332.5 keV are 
emitted in more than 99.9% of decays giving 2506 keV 
of deposited energy when both gammas get fully ab- 
sorbed in the NC volume. Once the PMT parame- 
ters and the discriminator threshold were found with 
^^Co, a fine tuning of the gain was done with the cen- 
tral LED. Then, an AmBe source was used to verify 
the ability of the system to detect neutrons. An ad- 
vantage of this neutron source is that in 60% of the 
cases the (a,n) reaction on Be leads to an excited 
^^C* state, which is then transfered to the ground 
state within 61 fs by emitting a gamma of 4438 keV. 
This gamma is a perfect trigger to start the acquisi- 
tion window in order to look for a following neutron 
capture signal. In this respect it nicely resembles the 
primary trigger during a physics run and gives an 
additional reference for the energy calibration. It is 
also to be expected that some of the neutrons from 
AmBe source will first produce a prompt recoil pro- 
ton and only then be captured on Gd. Thanks to 
the fiash ADC used one can separately look into an 
energy spectrum of the first hit and any following 
one within the same event. As it can be seen from 
Fig. [7j the primary hits form predominantly a bell- 
shaped spectrum specific to a prompt gamma while 
the secondary hits lead to a smoother spectrum which 
extends to higher energies, as one expects for par- 
tially escaping gammas appearing due to the neu- 



tron capture (see Sect. 3.2). The same behaviour is 
also confirmed in GEANT4 simulations (toolkit ver- 
sion 9.2.p01) of AmBe source where elastic scatter- 
ing off hydrogen, carbon and gadolinium nuclei, in- 
elastic scattering on carbon, capture on gadolinium 
and hydrogen nuclei and an energy quenching present 
for less ionizing particles are taken into account 
(FigjTj crosses). The initial kinetic energy of neu- 
trons is generated according to an energy distribution 
of AmBe source. A simplified model with a common 
fixed threshold and an energy resolution scaled as 
AE{E) = AE{A.4 MeV) x y'E (MeV) /4.4 (MeV) 
is used in the simulations, with a 2.8 MeV threshold 
and the energy resolution A£'(4.4MeV) = 700 keV as 
best fit. While the same values of the threshold and 
the resolution are used in the simulations to fit both 
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Table 1: Capture time simulated in GEANT4 as a function 
of Gd content. In order to extract Tcapture the same interval 
and binning as in the measured data were used for fitting. No 
background is assumed. 



Gd content, % 


0.18 


0.19 


0.20 


0.21 


0.22 


'^capture •> 


17800 


17100 


16290 


15600 


15000 


^'Tcapture-, 


±100 


±100 


±90 


±80 


±80 



measured ADC hits, the energy scaling of two spectra 
differs by 7%. This is due to a position dependence of 
the light collection and the fact that 4.4 MeV 7's are 
absorbed closer to the place of the AmBe source than 
neutron-capture gammas. Although the light propa- 
gation is not yet taken into account, a general good 
agreement between the measurement and the sim- 
ulations is obtained (FigjTj). Another verification of 
neutron capture on Gd is possible with the timing in- 
formation, i.e. by looking into time intervals between 
the start of the event and any following hit (Fig. [s]). 
Neutrons of several MeV need first to thermalize in 
the medium of the detector and only then get cap- 
tured following an exponential law with a characteris- 
tic time, the capture time, depending on the Gd con- 
tent of the used scintillator. Measured data (Fig. [s]) 
confirm well this description. Fitting the falling part 
of the spectrum and considering a flat background, 
one extracts Tcapture = 16.6 ± 0.3 /is. The same time 
constant is obtained in GEANT4 (Tab.[T]) where var- 
ious values of the Gd content can be considered. The 
simulated data are in good agreement with the mea- 
sured Tcapture and 0.20% of Gd admixture specified by 
the manufacturer (Tab. [T]). For further comparison 
we use the mean time of the measured distribution, 
i^) AmBe ~ ^^-^ A's, bccause it accounts also for the 
effect of neutron moderation and can be easily used 
in case of low statistics. 

The calibration runs using the ^"^Co and AmBe 
sources are regularly executed in order to monitor 
the response of the detector under well-defined cir- 
cumstances. 




TDC #ch (ns) 

Figure 8: AmBe measurement: time difference between a start 
of the event (first hit) and any following TDC hit within the 
acquisition time. The falling edge is fit with an exponential 
function and a flat background component. 

5 Preliminary physics data 

As a proof-of-principle we show in the following a pre- 
liminary analysis of physics runs. The muon-induced 
neutrons can be identified by the following features 
as described before: they produce multiple correlated 
hits within one event with distribution of time inter- 
vals similar to the one obtained for the AmBe source 
(Fig. [s]). The energy distribution of secondary ADC 
hits extends to higher values, i.e. well above the back- 
ground spectrum. In the first analysis, we identify 
a prompt muon together with correlated secondary 
neutrons as to study the neutron production by the 
well-tagged muons. The muon candidates are defined 
by a coincidence of the //Veto module on top of the 
neutron detector and the muon PMT's. Then, by 
using the time stamps of these events, a correspond- 
ing set of NC data is selected. Fig. [9]^Ze/i) shows the 
measured number of secondaries following a muon 
event and an expectation of Poisson distributed back- 
ground of accidental coincidences based on the single 
hit rate. Only those secondaries are accepted which 
are at least 200 ns apart from the primary hit and 
each other, in order to avoid an effect of PMT after- 
pulses. As can be seen, there is a clear excess of mea- 
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Figure 9: First data from a physics run, events selected with the /iVeto module on top of the NC and the muon PMT's first: 
number of secondaries within one event vs. Poisson expectation (left) and the time intervals between first and any following 
TDC hit (right). Live time of both systems running together is 51 days. 
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Figure 10: Preliminary data from a physics run, only neutron counter selected events: number of secondaries within one event 
vs. Poisson expectation (left) and the time intervals between first and any following TDC hit (right). Live time of the NC for 
this data set is 56 days. 
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sured secondaries above the background (A-values on 
Fig. |9]^/e/i)). The corresponding time distribution of 
these secondary hits (Fig.^right)) has a mean value 
(t) = 17 lb 3 /xs, which is clearly different from a 
central value of 29.8 /is specific for a random hit dis- 
tribution within the DAQ window of 59.6 /US, and in 
agreement with {t) j^^^^ = 18.5 /is. Furthermore, we 
performed a Kolmogorov-Smirnov (KS) test on un- 
binned data of Fig. [sjand Fig. [9]^ng/ii). This yields 
prob.(KS)=0.1 showing consistency of both spectra 
at a significance level of a=0.05. In overall, within 
the 51 live days analysed here we identify 24 neu- 
trons detected after the muons with an accidental 
background expectation of 0.2. This gives us a count 
rate of about 0.5 muon-induced neutrons per day. 

To search for neutrons produced due to untagged 
muons, we look into neutron counter data only and 
multiple hit events. This is possible if the neutron 
PMT's are tuned in such a way that the rate of acci- 
dental coincidences is well below an expected rate 
of muon-induced neutrons. To fulfill this require- 
ment, half of the neutron PMT's was set to a re- 
duced gain. Fig. ITO] presents corresponding distri- 



butions of the number of secondaries and of related 
time intervals. As in the first analysis, the time cut 
of 200 ns is applied between consecutive hits. As one 
can see, there is again a clear excess of measured data 
above the Poisson-based expectation of the back- 
ground (A-values on Fig. |10|^/e/;t)). The mean value 
(t) = 18.7 lb 1.7 fis is also in good agreement with 
i^)AmBe ~ ^^'^ t^^' above, KS-test is applied on 
unbinned data of Fig. [s] and Fig. [To| rig/i^) and re- 
sults in prob.(KS)=0.3, denoting consistent spectral 
shapes at a significance level a=0.05. The live time 
of this data set is 56 days. The data selection re- 
sults in 75 events with minimum one secondary pro- 
duced, among which there are 10 events with 2 or 
more secondaries detected. The expected acciden- 
tal background is 7.8 and contributes to events with 
only one secondary particle. First simulations show 
that correlated background from ambient neutrons 
contributes less than 10 events and from the Bi-Po 
decay chain less than 1 event, both with 1 secondary 
only. The full sample was also checked for a coinci- 



dence with the /xVeto module on top of the neutron 
detector and the muon PMT's. Ten events with at 
least one secondary could be attributed to a tagged 
muon. Thus, the excess shown in Fig. |10[/e/;t) clearly 
indicates muon-induced neutrons without requiring 
to detect the primary muon. 

It is evident, that the first physics data show qual- 
itatively the ability of the system to register muon- 
induced neutrons while a quantitative analysis of the 
production rate of these neutrons requires longer time 
of running to gather larger statistics and further de- 
tailed simulations to extract accurate detection effi- 
ciencies. 



6 Conclusion and outlook 



A neutron detection system based on 1 m^ of Gd- 
loaded scintillator has been described. The purpose 
of this dedicated instrument is to measure specifically 
the muon-induced neutron background in the envi- 
ronment of present and planned dark matter experi- 
ments, e.g. EDELWEISS-II and EURECA, hosted 
by the LSM underground laboratory. This back- 
ground is of concern for the DM searches and thus 
the neutron detector is in-location of the DM setups. 
The NC system was installed in fall 2008 and, as 
the calibration and first physics runs show, it is fully 
operational. The aim of the measurement is to de- 
tect, in particular, the neutron yield in Pb both with 
or without tagging the primary muons. More de- 
tailed and accurate Monte-Carlo simulations taking 
into account all the aspects, e.g. the muon interac- 
tion, the neutron production, a light propagation, are 
under further development. The preliminary analy- 
sis of physics data show that the count rate of muon- 
induced neutrons in the order of 1 neutron/day has 
been achieved. Data-taking will continue through 
2010 and early 2011 to get sufficiently large statis- 
tics. 
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